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Abstract 

Starting from the non-relativistic Pauli description of spin i particles, a set of fluid 
equations, governing the dynamics of such particles interacting with external fields and 
other particles, is derived. The equations describe electrons, positrons, holes, and similar 
conglomerates. In the case of electrons, the magnetohydrodynamic limit of an electron-ion 
plasma is investigated. The results should be of interest and relevance both to laboratory 
and astrophysical plasmas. 

PACS numbers: 52.27.-h, 52.27.Gr, 67.57.Lm 

1 Introduction 

The concept of a magnetoplasma has attracted interest ever since first introduced by Alfven [1], 
who showed the existence of waves in magnetized plasmas. Since then, magnetohydrodynamics 
(MHD) has grown into a vast and mature field of science, with applications ranging from solar 
physics and astrophysical dynamos, to fusion plasmas and dusty laboratory plasmas. 

Meanwhile, a growing interest in what is known as quantum plasmas has appeared (see, 
e.g., [2,3]). Here a main line of research can be found starting from the Schrodinger description 
of the electron. Assuming that the wave function can be factorized 0, one may derive a set of 
fluid equation for the electrons, starting either from an A^-body description, a density matrix 
description, or a Madelung (or Bohm) description of the wave function(s) [2,4]. As in classical 
fluid mechanics, the set of equationsn may be closed by a suitable assumption concerning the 
thermodynamical relation between quantities. These descriptions of the electron fluid, and 
its interaction with ions and charged dust particles, has been shown to find applications in 
many different settings [5-14]. Bart of the literature has been motivated by recent high field 
experimental progress and techniques [15-17], or the emergence of new areas, such as spintronics 
[18]. 

Indeed, from the experimental perspective, a certain interest has been directed towards the 
relation of spin properties to the classical theory of motion (see, e.g., [19-31]). In particular, 
the effects of strong fields on single particles with spin has attracted experimental interest in 
the laser community [21-26] . However, the main objective of these studies was single particle 
dynamics, relevant for dilute laboratory systems, whereas our focus will be on collective effects. 
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Moreover, strong external magnetic fields can be found in astrophysical environments such as 
pulsar [32,33] and magnetars [34]. Therefore, a great deal of interest has been directed towards 
finding good descriptions of quantum plasmas in such environments [35-38]. Thus, there is 
ample need and interest in developing models that are suitable for a wide range of applications, 
taking into account collective effects in multi-particle systems. There are in the literature several 
different studies of the effect of spin in plasma systems, e.g. kinetic descriptions in connection 
to fusion studies [39,40] and electron spin waves in solid-state plasmas [41]. Moreover, the 
connection between microscopic and macroscopic spin dynamics in the plasma state has been 
discussed by de Groot and Suttorp [42] . 

Inspired by both the historic and recent progress on quantum plasmas, a complete set of 
multi-fluid spin plasma equations was presented in Ref. [3]. In the current paper, we show, 
starting from the non-relativistic Pauli equation for spin ^ particles, how a set of plasma equa- 
tions can be derived for such spin ^ particles. These particles may constitute electrons, positrons 
(albeit non-relativistic), holes, or similar. Allowing these to interact with ions or charged dust 
particles, as well as other spin ^ particles, gives the desired governing dynamics of spin plasmas. 
We furthermore derive the appropriate magnetohydro dynamic description for such quantum 
plasmas, and investigate the effects of spin on the dynamics of the plasma. The limitations and 
suitable parameter ranges of the derived governing equations are discussed. The results should 
be of interest for both laboratory and astrophysical plasmas. 



2 Governing equations 

In relativistic quantum mechanics, the spin of the electron (and positron) is rigorously introduced 
through the Dirac Hamiltonian 



H = ca ■ (p+ -A^j - e4> + (5m e c 2 , (1) 



where a = (01,02,03), e is the magnitude of the electron charge, c is the speed of light, A is 
the vector potential, <f> is the electrostatic potential, and the relevant matrices are given by 

Here I is the unit 2x2 matrix and <x = (cxx, 03, 03), where we have the Pauli spin matrices 

ai = ( 1 ) ' a2 = ( i ~o ) ' and a * = ( I -1 ) • (3) 

From the Hamiltonian (pQ), a nonrelativistic counterpart may be obtained, taking the form 

Thus, the electron possesses a magnetic moment m = — /j,B{t(;\(T\il>) / , where [is = eh/2m e c 
is the Bohr magneton, giving a contribution — B ■ m to the energy. The latter shows the 
paramagnetic property of the electron, where the spin vector is anti-parallel to the magnetic 
field in order to minimize the energy of the magnetized system. According to (Jl|) and the 
relation dF/dt = dF/ dt + (1/ih) [F, H] , where F is some operator and [, ] is the Poisson bracket, 
we have the following evolution equations for the position and momentum in the Heisenberg 
picture [43,44] 

^ dF = ^( p+ c A )< (5) 



2 



dt 

while the spin evolution is given by 



where the spin operator is given by 



m e ^ = -e(E + -xB)-^^ B V(B-S), (6) 



It = H» bB X * (7) 



r- (8) 



The above equations thus gives the quantum operator equivalents of the equations of motion for 
a classical particle, including the evolution of the spin in a magnetic field. 

Next, we derive a set of nonlinear fluid equations for a gas of interacting electrons/positrons/holes 
(for a discussion of a kinetic approach, see Refs. [39] and [40]). The non-relativistic evolution of 
spin | particles, as described by the two-component spinor is given by (see, e.g. [4]) 



ih- 



(a) 



dt 



n 2 



2m j{a) 



#(«) (9) 



where mi a -\ is the particle mass, A is the vector potential, the particle charge, /X( Q ) the 
particle's magnetic moment, a = (<ti, a 2, 03) the Pauli spin matrices, 4> the electrostatic poten- 
tial, and (a) enumerates the wave functions. For an electron, the magnetic moment is given by 
\x e = —fJ,B = —eh/2m e c. From now on, we will define fj,r a \ = \x = qh/2mc. 
Next we introduce the decomposition of the spinors as 

&(a) = y^Ha) exp(i5( a )/fi)^( a ), (10) 

where n( a ) is the density, S7 a ) is the phase, and ipt a \ is the 2-spinor through which the spin ^ 
properties are mediated. Multiplying the Pauli equation ([9]) by , inserting the decomposition 
(|10p . and taking the gradient of the resulting phase evolution equation, we obtain the continuity 
and moment conservation equation 



dn {a) 
dt 

and 



+ V- {n (a) v (a) ) =0 (11) 



??? 



(a) + V H ■ VJ v {a) = q {a) {E + v (a) x B) 
+^(V ® B) • a (a) - VQ (a) V • (n (a) Z (a) ) (12) 

respectively. The spin contribution to Eq. (|12p is consistent with the results of Ref. [42]. Here 
we have introduced the tensor index a, b, . . . = 1, 2, 3, the velocity is defined by 

v (a) = — (vS (a) - ihip^Vcp) - —^—A, (13) 
m {a) V J m {a) c 

the Schrodinger like quantum potential (or Bohm potential) is given by 



n 2 



Q(*) = ~ - - i/2 v2n liy ( 14 ) 
2m W"(«) 
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the spin density vector is 

s (a) = ^l^fia), (15) 

which satisfies \st a \ \ = h/2, and we have defined the symmetric gradient spin tensor 

*(„) = (Vs(a)a) ® (Vflf a) ). (16) 
Moreover, contracting Eq. by •P , Jn<t, we obtain the spin evolution equation 

+ «(«) ' a (a) = -^B x + — -L— S(a) x [d a (n (a) d a s {c0 )] . (17) 

We note that the particles are coupled via Maxwell's equations. 

Suppose that we have N wave functions for the same particle species with mass m, magnetic 
moment fj,, and charge q, and that the total system wave function can be described by the 
factorization \P = !Ph)<P(2) • • • *(2V)- Then we define the total particle density for the species with 
charge q according to 

N 
(a)=l 

where p a is the probability related to the wave function $? a ). Using the ensemble average 
(/) = YliaPa( n (a)/ n q)f (f° r an y tensorial quantity /), the total fluid velocity for charges q is 
V q = (vt a -\) and the total spin density is S = (s( a -\). From these definitions we can define the 
microscopic velocity in the fluid rest frame according to wr a -\ = vr a \ — V, satisfying = 0, 

and the microscopic spin density Sr a \ = sr a \ — S, such that (5( a )) = 0. 
Taking the ensemble average of Eqs. (fTTj) . (fT2j) . and (fTTj) . we obtain 



dn, 
~~di 



q + V ■ (n q V q ) = 0, (19) 



and 



mn i (jf t + y i- V ) V i = qn q (E + V q x B)-V -n q - VP q + C qi + F Q (20) 
% (§- t + V q ■ V) S = -^B xS-V-K q + n s (21) 



respectively. Here we have added the collisions C q i between charges q and the ions i, denoted 
the total quantum force density by 



^*(V ® B) ■ S - n q (VQ {a) ) - - V • KI ) - - V • (n g E ' 

~V ■ [n q (VS a ) <(VS ( a a) )> +n q ((VS {a)a )} (VS a )], (22) 

consistent with the results in Ref. [42], and defined the the nonlinear spin fluid correction 
according to 

fl s = -Sx [d a (n q d a S)} + -Sx [d a (n q (d a S {a) ))} 
to m v ' 

^ /fW x [da{ \ d a S)] + ^ /«W x [S a („ (a) fl- 5(a) )]\ , (23) 
to \n (a ) 7 m \ n (a) / 
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where n = mn[(w^ <g> w^) — \(w^)/3] is the trace- free anisotropic pressure tensor (I is the 
unit tensor), P = mn(u>? n) is the isotropic scalar pressure, Z = (V5" a ) (g) (VS* a ) is the nonlinear 
spin correction to the classical momentum equation, Z = ((V<S( a ) a ) <8> (ViS^)} is a pressure 
like spin term (which may be decomposed into trace-free part and trace), K = n(tO( a ) <S> S( a )) 
is the thermal-spin coupling, and [(V <g> B) ■ S] a = (d a Bi,)S b . Here the indices a, b, . . . = 1, 2, 3 
denotes the Cartesian components of the corresponding tensor. We note that the momentum 
conservation equation ([20]) and the spin evolution equation (}2~Tj) still contains the explicit sum 
over the N states. 

The coupling between the quantum plasma species is mediated by the electromagnetic field. 
By definition, we let H = B/po — M where M = 2n q pS/h is the magnetization due to the 
spin sources. Ampere's law V X H = j + to^tE then takes the form 

1 dE 

VxB = hqU + hi) + ^2^, (24) 

where we have the magnetization spin current jm = V x M and the free current j . The system 
is closed by Faraday's law 

Vx E =-_. (25) 
3 Electron— ion plasma and the magnetohydrodynamic limit 

The preceding analysis applies equally well to electrons as holes or similar condensations. We 
will now assume that the quantum particles are electrons, thus q = — e, where e is the magnitude 
of the electron charge. By the inclusion if the ion species, which are assumed to be described 
by the classical equations and have charge Ze, we may derive a set of one- fluid equations. The 
ion equations read 

_i + V ■ irnVi) = 0, (26) 

and 

d_ 

dt 



mini ( — + Vi ■ V ) Vi = Zen, (E + V t x B) - V • fl, - VP, + C iq . (27) 



Next we define the total mass density p = {m e n e + mini), the centre-of-mass fluid flow veloc- 
ity V = (m e n e V e + miniVi) / p, and the current density j = —en e V e + ZeniVi. Using these 
denfinitions, we immediately obtain 

^ + V • ( P V) = 0, (28) 

from Eqs. (I19p and (I26p . Assuming quasi-neutrality, i.e. n e ~ Zni, the momentum conservation 
equations (f20j) and ([27]) give 

p(j t + V -V^jV = j xB-V -U-VP + Fq, (29) 

where n is the tracefree pressure tensor in the centre-of-mass frame, P is the scalar pressure in the 
centre-of-mass frame, and the collisional contributions cancel due to momentum conservation. 
We also note that due to quasi-neutrality, we have n e = pj (m e + nn/Z) and V e = V — niij/Zep, 
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and we can thus express the quantum terms in terms of the total mass density p, the centre-of- 
mass fluid velocity V, and the current j. With this, the spin transport equation (|2ip reads 

| + V.v) S ^VS-^ B *^ + f)v.K fl+ (,„ c + f 



In the momentum equation (|29p . we have the force density j x B + Fq, where j is the current 
produced by the free charges. In general, for a magnetized medium with magnetization density 
M, Ampere's law gives the free current in a finite volume V according to 

j = — V x B - V x M, (31) 

where we have neglected the displacement current is jjj = eodtE. The delta-function contribu- 
tion of the surface current is an important part of the total current when we are interested in 
the forces on a finite volume, as will now be shown. 

It it worth noting that the expression of the force density in the momentum conservation 
equation can, to lowest order in the spin, be derived on general macroscopic grounds. Formally, 
the total force density on a volume element V is defined as F = limy_>o(X) Q fa/V), where f a 
are the different forces acting on the volume element, and might include surface forces as well. 
For magnetized matter, the total force on an element of volume V is then 



/tot = f jtot x BdV + / (M xnjxBdS 

JV JdV 



(32) 



where (neglecting the displacement current) j tot = j + V x M. Inserting the expression for the 
total current into the volume integral and using the divergence theorem on the surface integral, 
we obtain the force density 

F tot =jxB + M k VB k , (33) 

identical to the lowest order description from the Pauli equation (see Eq. (|29p ). Inserting the 
free current expression (f3Tj) . due to Ampere's law, we can write the total force density according 
to 2 

F i = -d l ( — -M B ) +d k (H i B k ). (34) 



, 2 Mo 

The first gradient term in Eq. (|34p can be interpreted as the force due to a potential (the energy 
of the magnetic field and the magnetization vector in that field), while the second divergence 
term is the anisotropic magnetic pressure effect. Noting that the spatial part of the stress tensor 
takes the form T ik = -H i B k + (B 2 /2fi - M ■ B)5 ik [42], we see that the total force density 
on the magnetized fluid element can be written F l = —dkT lk , as expected. Thus, the Pauli 
theory results in the same type of conservation laws as the macroscopic theory. The momentum 
conservation equation (|29p then reads 

K^ + F ' V ) F = " V (C- M ' B ) +i? ^- VP ' (35) 
where for the sake of clarity we have assumed an isotropic pressure, dropped the displacement 
current term in accordance with the nonrelativistic assumption, and neglected the Bohm poten- 
tial (these terms can of course simply be added to (i35|) ). 

Approximating the quantum corrections, using L Xp where L is the typical fluid length 
scale and Xf is the Fermi wavelength for the particles, according to [2] 

(VQ (Q) ) « -V ( — ^— VX /2 J = VQ. (36) 
\2mn q J 
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We then note that even if Q is small, the magnetic field may, through the dynamo equation 
(|25p . still be driven by pure quantum effects through the spin. 

A generalized Ohm's law may be derived assuming C e i = en e r/j, where r\ is the resistivity. 
From the electron mometum conservation equation (|20p combined with Faraday's law we obtain 



dB 

~dt 



V x B 



jxB 



en,, 



rjj + 



m e dV e 
e dt 



Fq 

en e 



(37) 



and 



where n is the resistivity. Here we have omitted the anisotropic part of the pressure, 
neglected terms of order m e /ra, compared with unity. 

The electron inertia term is negligible unless the electron velocity is much larger than the ion 
velocity. Thus whenever electron inertia is important, we include only the electron contribution 
to the current, and use Ampere's law to substitute V e = V xB/en e /iQ into the term proportional 
to m e in (|37h . which gives the final form of the Generalized Ohm's law 



dB 

~dt 



V x <V x B 



jxB 



en,, 



VJ 



e 2 /x 



0_ 

dt 



V x B 




In the standard MHD regime the Hall term and the electron inertia term are negligible. 
During such conditions the quantum force is also negligible in Ohm's law, which reduces to its 
standard MHD form 

_ = V x (V x B - rjj) (39) 

Note, however, that the quantum force including the spin effects still should be kept in the 
momentum equation (|29p . Equations (|28p . (|29p. and (|38p together with the spin evolution 
equation (i30l) . which is needed to determine Fq, constitutes the basic equations. In order to 
close the system, equations of state for the pressure as well as for the spin state are needed, as 
will be discussed in section 5. 



4 Closing the system 

The momentum equation, Ohm's (generalized) law and the continuity equation need to be 
completed by an equation of state for the pressure. As is well-known, rigorous derivations of 
the equation of state is only applicable in special cases of limited applicability to real plasmas. 
Useful models include a scalar pressure where the pressure is proportional to some power of the 
density, i.e. 

(4o) 

where d/dt = d/dt + V, • V is the convective derivative and 7 S is the ratio of specific heats, 
which in general can be different for different species s. Secondly, we note that the magnitude 
of the terms that are quadratic in the spin depends highly on the spatial scale of the variations. 
In MHD, the scale lengths are typically longer or equal to the Larmor radius of the heavier 
particles, which means that the terms that are quadratic in S can be neglected in the expression 
for the quantum force Fq as well as in the spin evolution equation (I30p . To lowest order, the 
spin inertia can be neglected for frequencies well below the electron cyclotron frequency. Also 
omitting the spin-thermal coupling term, which is small for the same reasons as stated above, 
the spin-vector is determined from 

B xS = 0. (41) 
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which has a solution 

consistent with standard theories of paramagnetism, as the spin anti-parallel to the magnetic 
field minimizes the magnetic moment energy, so that 

M = 2^ B n eV (^^jB. (43) 

Here B denotes the magnitude of the magnetic field, r/(x) = tanhx is the Brillouin function!! 
and B is a unit vector in the direction of the magnetic field. In this approximation the spin 
evolution equation (|30p can be dropped, and the quantum force can be written 



Q 



where the second term comes from the spin. Combining the approximations presented in this 
section together with the MHD equations presented in section 3, or either of the two dust systems 
presented in section 4, closed systems are obtained. 



5 Illustrative example 

Let us next consider the spin- model described by (28), (35), and (39) with the resistivity put 
to zero. Furthermore, we note that for most plasmas [ibB <C k B T, and thus we can use 
the approximation t](iibB /k B T) ~ fJ, B B /k B T '■ For definiteness, we also choose an isothermal 
equation of state such that VP = k B TVn. Next we let B = B Q z + B%, M = M z + Mi, 
n = no + ni, where index denotes the equilibrium part and index 1 denote the perturbation, 
and we have assumed that the equilbrium part of the velocity is zero. Linearizing around the 
equilibrium and Fourier analysing, and omitting the non-spin part of the quantum forcejl we 
find that the general dispersion relation can be written 



UJ 



k 2 C\ 



k 2 C\ + k^c 2 . 



- + Kk 







(45) 



where c s = [(ksTe + ksTi) /mi] 1 / 2 is the standard acoustic velocity. Formally, Eq. (|45f) looks 
similar to the standard ideal MHD dispersion relation. However, we note that firstly spin effects 
appear in the spin-modified Alfven velocity 



C A = C A [1 



mc 2 k B T 



1/2 



(46) 



1 /2 

where to pe is the electron plasma frequency and Ca = {Bq/hqPo) is the standard Alfven 
velocity. Moreover, another spin- modification is introduced in the factor c s , given by 



k B T 



1/2 



(47) 



2 In general, the thermodynamic equlibrium distribution of spin results in a magnetization proportional to the 
Brillouin function Bj(fiB /ksT), where j is the spin of the particles in question. For the special case of spin 
1/2-particles we have S 1 / 2 (a;) = tanh(s). 

3 The ordinary part of the quantum force is negligible compared to the spin coupling provided eBo ^> hk 2 . 
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The first factor of Eq. (|45p describes the spin-modified shear Alfven waves, and the second factor 
describes fast and slow magnetosonic waves. We point out that the assumption made in this 
section initially, fisB <C ksT, implies that both Ca and c s given by (j4"6j) and ([17J) respectively 
are real. Moreover, we stress that we cannot simply obtain the spin-modified magnetosonic 
dispersion relation from the classical one by the substitutions Ca — > Ca an d c s — > c s , since the 
classical acoustic velocity c s still appears in the factor (lo 2 — k 2 c 2 ), as seen from (|45p . We end 
this section by noting that the modification of the Alfven velocity can be appreciable for a high 
density low temperature plasma. 

6 Summary and Discussion 

In the present paper we have derived one-Fluid MHD equations for a number of different plas- 
mas including the effects of the electron spin, starting from the Pauli equation for the individual 
particles. In particular we have derived spin-MHD equations for an electron-ion plasma. Fur- 
thermore, the general equations derived in section 2 constitutes a basis for an electron-positron 
plasma description including spin effects. 

In order to obtain closure of the system, our equations needs to be supplemented by equations 
of state for the pressure as well as for the spin pressure. In the MHD regime, a rather simple way 
to achieve this closure has been discussed in Section 4, where we assume that the scale lengths 
are long enough such terms that are quadratic in the spin vector as well as the spin-thermal 
coupling are neglected. However, we here note that if more elaborate models are used, the spin 
pressure together with the spin-thermal coupling might play an important role in the generalized 
Ohms's law fl38l) . 

Since the spin-coupling give raise to a parallel force (to the magnetic field) in the momentum 
equation, the parallel electric field will not be completely shielded even for zero temperature, 
contrary to ordinary MHD. As an immediate consequence, the spin-coupling can give rise to a 
rich variety of physical effects. In this paper we have limited ourselves to present a single example, 
linear wave propagation in an electron-ion plasma, and shown the modification of the dispersion 
relation due to the spin effects. In particular, we note that the spin effects are important for 
low temperature, high densities, and/or strongly magnetized plasmas. The latter can be found 
in astrophysical systems, such as pulsars or magnetars [38], where however other modifications, 
such as vacuum polarization and magnetization [15,45], of our set of governing equations may 
be necessary in order to accurately model such plasmas. Moreover, ultra-cold plasmas may 
even be found in laboratory environments, where currently mK temperature plasmas can be 
formed [46]. Studies involving the spin dynamics through the spin evolution equation (12ip . 
kinetic effects associated with the spin, as well as nonlinear spin dynamics are projects for 
future work. Finally, we note that dusty plasmas can sustain weakly damped modes with low 
phase velocities [47], and quantum and spin effects tend to be important in this regime. 
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